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SUMMARY 

A study of the statistics of the look-to-look variations in the returns re- 
corded in-flight by a digital* realtime SAR system is being conducted for the 
Avionics Laboratory, Radar Branch, WPAFB, under USAF Contract F33615-77-C 
1169* The study objectives include the determination of whether these varia- 
tions demonstrate class-unique information content adequate to allow the use 
of some measure of that variance in discriminating between or identifying any 
or all of various classes of objects and features mapped by the SAR system. 
This paper briefly discusses the data, their source, the techniques employed 
in sampling, and some measures under study and presents a demonstration of 
the positive results thus far obtained toward meeting the objective defined 
above. The determination that the variations in the look-to-look returns, 
from different classes do carry information content unique to the classes is 
illustrated by means of a model based on four variants derived from the four 
look in-flight SAR data under study. The model was limited to four classes 
of returns: mowed grass on a athletic field, rough unmowed grass and weeds 

on a large vacant field, young fruit trees in a large orchard, and metal 
mobile homes and storage buildings in a large mobile home park in Reedley, 
California. A data population in excess of 1000 returns, distributed between 
the four classes, was selected. These returns represent over 250 indi- 
vidual pixels from the four classes. The multivariant discriminant model 
operated on the set of returns for each pixel and assigned that pixel to one 
of the four classes, based on the target variants and the probability dis- 
tribution function of the four variants for each class. The model, the 
variants used, and the test results are presented. 
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1.0 OBJECTIVES AND GENERAL APPROACH 


1.1 OBJECTIVES 

The overall objective of the study is to analyze the look-to-look variations 
of the multilook SAR returns from selected objects and features to determine 
whether some measure of that variance can be used in discriminating between 
or identifying any or all of the features and objects analyzed. 

1.2 GENERAL APPROACH 

The general approach is to analyze the variations of the returns from each of 
a number of sets of actual multilook returns representing statistically 
significant samples of selected classes and subclasses of known objects and 
features . 

These sample sets are comprised of the radar returns associated with single 
pixels identifiable with known objects and features. The single pixels are 
selected interior to gross features, scene segments, or objects (i.e., indi- 
vidual trees within an orchard or grove, individual structures within a 
group or within a cluster of structures or trees, etc. ) 

Initial designation of classes and subclasses to be sampled was based on a 
requirement to obtain a realistic spread of variations and power levels. 
Additional classes and subclasses have been added as the study has progressed. 
All selections are based on ground truth information to minimize the number 
of "wild points" which may occur. 

2.0 DATA SOURCE AND THE SELECTION PROCESS 

2.1 SOURCE OF DATA 

The SAR returns used in the study are obtained from digital data recorded 
in flight on wideband magnetic tapes during the recently concluded, USAF- 
conducted, 72-flight experimental flight test series with the in-flight, real- 
time data processing FLAMR ( Forward Looking Airborne Multimode Radar ) sys- 
tem . This system operated in the 1 6 to lb. 5 GHz band . 

Approximately 1+50 wideband, 30.5 x 2.5 cm ( lh x 1 in.), ll+-track magnetic 
tapes were recorded during the flight test program. These tapes and asso- 
ciated data and equipment are now located in the SAR Data Bank at the Applied 
Research Laboratories, The University of Texas at Austin (ARL:UT). The Data 
Bank is maintained for purposes of supplying SAR data to qualified. Air Force- 
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approved members of the technical and scientific community. The Data Base is 
supported by the Air Force through USAF Contract F33615-77-C-1169 with ARL:UT. 

2.2 THE SELE CTION PROCESS 

The flight data to be examined are identified by review of the documentation, 
the available ground truth data* and photographs of the imagery associated 
with various flights. 

The next step is to replay the appropriate wideband tape corresponding 
to the selected leg of a flight. This replay is accomplished by means of the 
wideband tape drive shown on the right in the photograph in Fig. 1. The 
recorded imagery is reviewed on the Monitor Display in 16 gray level format. 
The DBS imagery used in the study is shoim on the display in 360 vertical 
lines comprised of 38*4 range bins each. 
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FIGURE 1. GAR DATA BANK FLIGHT DATA RETRIEVAL EQUIPMENT 


When the desired map is painted on the display, the image is then "frozen" 
and the tape drive halted. The desired objects or features to be sampled are 
then located on the display. To locate the brightest pixel in the case of 
small discrete objects, the set of gray scale emphasis or deletion switches 
below the monitor is used. To obtain the azimuth line and range bin 
numbers of a selected pixel, the manual controls below the displa are 
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, used to position a local cursor over the pixel. The desired coordinates are 
then read directly from the cursor controls. 

After all desired sample pixels on the scan have been identified and coordi- 
nates noted, a CCT (Computer-Compatible Tape) is placed on the small tape 
drive visible at the left in Fig. 1. The data select /dump controls on the 
panel to the right of the display are then set to dump the desired data. 

The practice thus far in the study has been to dump the full scan of FM data 
(Filter Magnitudes recorded in flight at the output of the Doppler Processor) 

- and a file of the associated RIOT data (Radar Input-Output data recorded from 
the radar-computer I/O bus). The flight data cure in computer-compatible 
format and are processed on a general purpose computer. 

The FM data are next reformatted and five arrays of FM data, similar to those 
appearing in Fig. 2, are printed for each selected pixel for use in checking 
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FIGURE 2. FM PRINTOUT FOR 35-PIXEL SECTIONS OF THE FOUR 
SINGLE-FREQUENCY AND THE COMPOSITE 4:1 MAPS 

and finalizing the selection and coordinates for the sample pixels. The 
FM values appearing in Map 1234, the composite 4:1 overlaid map, are examined 
to verify that the selected coordinates of the center point of the printout 
array represent, for the case of a small discrete, the brightest return from 
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the discrete. Iij an immediately^ adjacent pixel has a higher FM value, then 
the coordinates are changed to the coordinates of the brighter pixel. Thus 
far in the study, no further effort to correct for range/Doppler straddling 
has been made. After the coordinates are finalized, they are used with com- 
puter programs and the two data tapes, FM and RIOT, to combine the returns 
data and other required information on a punched card for each sample pixel. 
The punched cards are then used as the data bake for analysis. 


2.3 THE TEST CASE 

A photograph of the imagery from which the sample sets treated in this paper 
were obtained is shown in Fig. 4. This imagery is one of the DBS scans 
made during a mapping flight in the San Joaquin Valley area of California 
near the town of Reedley. This scan was chosen for sampling because of the 
diversity of terrsin, vegetation, and cultural features and the large Data 
Bank stock of photographs of various portions of the town and nearby country- 
side. Figure 3(A' is a photograph of a moved grass athletic field located 
Just above and to the right of the two bridges visible in the upper left 
corner of Fig. 4. The samples comprising sample class 201 were obtained from 
this location on the imagery. 

Figure 3(B) is a photograph of a large vacant field covered with rough grass 
f ,and^Wiifcls and located to the right of the group of bright returns from the 
buildT&gs located Just left of center at the extreme top edge of Fig. 4. 

The samples from this field comprise class 202. It was anticipated that 
there would be substantial overlap between the variants for classes 201 and 
202, and the two similar features were selected for this reason. 

Class 102 is comprised of fruit trees from the orchard shown in Fig. 3(C). 
This orchard is the large triangular area located Just above the river in the 
lower right quadrant o^ Fig. 4. 

Class 9000 was selected from the rectangular clubcer of bright returns Just 
below the sharp bend in the river in the lower right corner of Fig. 4. These 
returns are from a large mobile home park, visible in Fig. 3(D). 
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Since the I and Q data available the Data Bank have not yet been studied, 
only the FM data will be discussed, here. As shown in the printout of FM 
data illustrated in Fig. 2, there are five arrays of 5 x T FM values each, a 
total of 35 pixels. Each array carries a MAP number and an FRQ value code to 
the right of the array, with the exception of the array at the bottom of the 
figure. The column of numbers on the left edge of the picture contains the 
range bin coordinate for the associated row of pixels, and the map line 
numbers above the top array denote the azimuth coordinate for the associated 
column of pixels. The arrays for maps 1, 2, 3, and 4 are obtained from the 
four complete 1:1 overlay ratio, single frequency maps formed in ■‘•’he process 
of producing the 4:1 overlay map from which the bottom array is obtained. 
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flote that each pixel with a given range bin and azimuth line set of 
coordinates is represented in all five arrays. 


ARMY No. 
LINE No. 
NAM 
FR£Q..f 0 


n- A rray- 1 Array-5 — * k- A rray-9 -h 

123456781234567812345678---- 

iimiiiiiiiiiiiiiiiiiii 


ARMY No. 
LINE No. 
HAM 
FREQ.-fj 


Array-2 -h Array-6 -H Array-10 -h 
123456781234567812345678 

minimimiiiiiiiiii 


ARMY No. 
LINE No. 
NAM 
FREQ.*f 2 


K- Array- 3 -H k- Array- 7 K- Array- 11 -h 

123456781234567812345678 

iiiiiiiiiiiiiiiiiiiiiii 1 


ARMY No. 
LINE No. 
HAM 

FRtQ.-f 3 


r Array-4 -h h- Array- 8 — h k— A rray- 12 — h 
1^3456781234567812345678 

1 1 ii 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


4:1 Overlaid 
LINE No. 
HAM 234 


1 2 3 4 5 6 7 8 9 10 12 14 16 18 20 22 24 

1 1 1 1 m 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


Rixa no. 

Return 

tlM sequence 


XX^XXXXXXXXXXXXXXXXXXX 

1 2 3 4 5 6 7 8 9 10 12 14 16 18 20 22 
f 0 f l f 2 f 3 f 0 f l f 2 f 3 f 0 f 1 f 2 

f l f 2 f 3 f 0 f 1 f 2 f 3 f 0 f l f 2 f 3 

f 2 f 3 f 0 f l f 2 f 3 f 0 f l f 2 f 3 f 0 

f 3 f 0 f l f 2 f • f 4 f l f 2 f 3 f 0 f l 


X X 

24 


FIGURE 5 - FORMAT I C ! U OF MAP OVERLAY LINES WITH 

ASSOCIATED FREQUENCY SEQUENCE OF RETURNS 
The FRQ value given for each of the 1:1 overlay maps represents the trans- 
mitter frequency use i in the formation of that map for this particular scan. 
This frequency varies f'rom flight to flight ar.d from scan to scan, but its 
cycle of occurrence does not vary. The cycle sequence is in order of descend- 
ing frequency (i.e., 16.U, l6.3, 16.2, l6.1, or any combination such as 16.3, 
l5.2, 16.1, 1 6.U, and so on). When the "frequency hep" mode was employed, as 
it often was, the transmitter frequency changed for each 8-azimuth-line 
array formed. 
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The Method of overlay employed in forming the separate 1:1 naps and the 
composite kil overlay nap is illustrated in Fig. 5- The key datum which the 
authors wish to convey here is the association of a particular frequency 
sequence with each pair of azimuth lines in the composite H:1 overlay map, 
which is farmed by averaging together the appropriate pixels from the four 
1:1 maps, 1, 2, 3, and b. Mote that three 8-line arrays are presented for 
each of the four 1:1 nape and that each array bears an array number, repre- 
senting the order in which the arrays are formed by the system. Array Ho. 1 
is formed first, then array Ho. 2 is formed but with an azimuth displacement 
of two azimuth lines to the right. Hext come arrays No. 3 and Ho. U, respec- 
tively, each with an additional 2-line displacement relative to the 
preceding array. 

After the first four arrays are formed, note that there are only two azimuth 
lines of the k:l map that can be formed by averaging together the corres- 
ponding lines of the four 1:1 maps. Thus, lines 7,8 of array Ho. 1, lires 
5,6 of array Ho. 2, lines 3, 1 * of array Ho. 3, and lines 1,2 of array Ho. U 
correspond in azimuth and, averaged together, form lines 1,2 of the U:1 map. 

The row of x's below the U:1 map denotes pixels on each map azimuth line. 

The column of f's below each pair of x's denotes the time order in which each 
return frequency will be received by the SAR. For pixels 1 and 2, the fifst 
transmitter frequency illuminating these points was f^, then, in order f^, 
fg, and f 2 * For the next pair of pixels the first frequency will be f^, 
from i.ines 7 and 8 of array No. 2. The rest of the sequence is f^ from 
lines 5,6 of array No. 3, fj from lines 3,1* of array No. U, and finally, f^ 
from array No. 5. These sequences continue to cycle throughout, the map. 

This varying frequency sequencing has been found to be quite important in the 
study of variants and the study of its effects is continuing. 

3.2 SCINTILLATION 

Four enlarged 1:1 map strips of the mobile home park and adjoining fields 
shown previously in Figs. 3(D) and U are presented in Fig. 6. Note the map- 
to-map variations in the pattern of returns from the structures in the park 
and the Surrounding fields. Gray level 8 was emphasized in these pictures 
and the higher gray shades were deleted to better illustrate the variations. 
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These variations are also apparent in the FM arrays printed in Fig. 2. 

Bote the variation in FM values between pixels with corresponding range 
bin and line coordinates. It has been discovered that the information con- 
tent in these variations is significant in feature discrimination as is 
demonstrated in the final sections of this paper. 

k.O SAMPLE VARIABTS 

Four variants were generated from the U-look FM data discussed previously. 

The four sequential values representing each pixel, converted to power units, 
were used to generate four variants employed in the classification model. 

The FM values are related to power by log^ P plus a constant so the 
conversion to power from the log filter magnitudes is given by 

P - K^ 8 ) , (1) 

where K is a calibration factor which may change from scan to scan or within 
a single scan. To illustrate the classification model, data were selected 
from regions within the scan over which K remained constant. 

The four target variants based on measurement are the mean, Z^; sample 
variance, Z^i fast variation, Z^; and slow variation, Z^, defined as 


z i - r ^ p i 


(2) 

<VV 2 


(3) 

Z 3 * [ I P 1- P 2 I * IV P 3 I 

+ |P 3 -P U I 1/3 

(U) 

■ tlv p 3 l * l p 2- p j 

+ lv p ul 1/3 

(5) 


It should be noted in the above relations that variants 1 and 2 do not 
involve time sequence of the data whereas the variants, fast variation and 
slow variation, do Involve the time sequence of the returns. Other measures 
of dispersion investigated, but not reported here include average deviation, 
color, deviation from color, total variation, spread, and slope. 
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5.0 STATISTICS OF THE VARIANTS 


5.1 DISTRIBUTIOHS INVESTIGATED FOR THE VARIANTS 

Since the probability distribution function of the measured data or variants 
derived from the measured data was required in the formulation of the classi- 
fication model, histograms of the distributions of the various variants were 
constructed from the sample data sets. Four standard probability distri- 
bution functions, which included the Lognormal, the Weibull, the Chi-Square, 
and the Rayleigh PDF's, were fitted to the data. It was found that the Log- 
normal and Weibull PDF's generally fitted the test data well with the Weibull 
POP giving a better fit than the Lognormal PDF. 

The Lognormal and the Weibull distributions are written as 

p(X) = — r expl-UnCx/M)) 2 ^^ (6) 

/2vXOjjj 

and 

p(X) = (X/o ) n_1 exp{-(X/o w rt (7) 

respectively, where M is the median of X, is the standard deviation of 
tn(X/M), o v is a scale factor, and n is a shape parameter. 

5-2 EXAMPLE HISTOGRAMS AND CURVE FITS 

Histograms of unnormalized power for samples of rough grass and mobile homes 
along with the corresponding Lognormal curve fits are shown in Fig. 7 for 
returns from individual SAR maps for each of the four looks with 1:1 overlay. 
It should be noted that there is little change in the distribution curves 
between maps for the two extremes of power levels represented by the two 
classes. The distribution and curve fits, both Lognormal and Weibull, are 
also shown in the same figure. The Lognormal, in general, appears to fit the 
histogram data but, for some data, does not appear to be veil centered about 
the histogram data. This apparent anomaly occurs whenever o TU is large (on 
the order of 1.50). Histograms and Lognormal distributions for received 
power for the four target types are shown in Fig. 8. These curves show that 
the Lognormal distribution is a good representation for the PDF for the four 
sample classes. 
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FIGURE 7 

RECEIVED POWER DISTRIBUTION - HISTOGRAMS AND FITTED LOGNORMAL FUNCTIONS 
FOR TWO TARGET TYPES SHOWING TYPICAL VARIATION BETWEEN MAPS. THE LOWER 
SET WHICH SHOWS THE 4-MAP MEANS ALSO SHOWS A COMPARISON BETWEEN THE WEIBULL 
AND LOGNORMAL PDF. 
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ATTERGRAMS OF FAST VARIATION IN POWER VERSUS MEAN POWER AND OF SLOW 
VARIATION IN POWER VERSUS MEAN POWER FOR FOUR CLASSES OF TARGETS 



To further illustrate the dispersion of the variants between sample classes, 
a sample set of crossplots of the variants, fast and slow variation, versus 
mean power is shown in Fig. 9 for the four classes.. 

6.0 A MULTIYABIABT DI3CRIMINAHT MODEL AHD RESULTS WHEN APPLIED TO A TEST 
CASE 

6.1 CHOICE OF THE CLASSIFICATION MODEL 

To demonstrate the utility of several of the variants currently being stud: ed 
for discrimination, a classification model based on the statistics of selected 
variants derived from the SAB return data was formulated and applied to the 
test sample. Four variants and four classes of objects were used to demon- 
strate the effectiveness of the Discriminant Model; however, it should be 
noted here that the model is by no means limited to these dimensions. 


Although the Weibull PDF (Probability Distribution Function) fitted the 
variants for all four classes comprising the test sample, the Shape Factors 
were large numbers for some targets, thereby making the Weibull PDF less 
suitable than the Lognormal PDF for construction of the classification model. 
Consequently, the Lognormal distribution was chosen for constructing the 
classification model. 

6.2 DESCRIPTION OF THE MODEL 

The variants, (Z^, Zg, Z^, Z^) , derived from on. . rt of measurements are 
assumed to belong to one of four sample classes comprising the test sam- 
ple. Since these quantities have a Lognormal distribution, the PDF of Z, 
when Z belongs to class , is represented by p( Z ) and for the vector Z is 
given by the following: 


P(Z|T.) 



( 8 ) 


where 


Z* [fcn(Z 1 /M n ) , tn(Zg/M i2 ) , ... ln(Z /M ±u ) ] 


Z' is the transpose of Z, and 


Z. *7T(Z /M, ) 
pi r ir 
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for r«l,2,3»4, are the medians of Z y when belongs to T^. Cj is the 
covariance matrix of Z. A decision function is formed that assigns each 
vector Z to one of the sample classes. This function will assign Z to class 
Tj if the following condition 

pUl^) > p(Z|Tj) 

is met for all J not including i. 

Inserting the expression 

p(z|T i ) 
r iJ " in p(Z|Tj) 

into the ratio of two probability density functions , one obtains a decision 
function as follows 


(9) 


( 10 ) 



A decision is now made, based on Eq. (11), as to which class Z belongs. By 
Eq. (9), Z is assigned to Class if r^>0, and is not assigned to if 


For the selected example given here, only the diagonal terms in the 

covariance "<at rices were retained. The off-diagonal terms were set to zero. 

All combinations for r^, were computed and Z was placed in class k when the 

condition r, ,>0 was met for all j. 

■KJ 

6.3 TEST RESULTS 

Employing the data representing the four sample classes described in 
section 4, the model was tested against each data point for several combina- 
tions of measurement vectors Z. The result of one test is shown in Table I 

when the four variants— Mean. , Standard Deviation,, . Fast Variation, and 

pvr pwr 

Slow Variation — were employed in the model. As is shown in Table I, data 
points of class 102 (young fruit trees in an orchard) were identified 
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TABLE I 


NUMBER OF CORRECT AND INCORRECT CLASSIFICATIONS 
FOR EACH PIXEL VERSUS EACH CLASS 


Class 

Class 

Description 

Number of 


Samples 

102 

201 

202 

9000 

102 

Fruit Trees 

50 

34 

5 

11 

0 

201 

Mowed Grass 
Field 

49 

0 

44 

5 

0 

202 

Unmoved Grass 
A Weed Field 

54 

1 

4o 

l4 

0 

■ 

Metal Mobile 
Homes 

37 

0 

0 

0 

37 


correctly by a ratio of slightly better than 2:1 over incorrect classifica- 
tions as belonging to the grass classes 201 and 202. The ratio of correct 
to incorrect classifications for the moved grass class (201) was nearly 9:1 
and here the incorrectly classified pixels were classified as belonging to 
the rough grass class. The ratio of correct to incorrect for the rough grass 
pixels (202) was about 1:4, a poor shoving. Hovever , the grasses vere 
initially selected to test the variants for tvo very similar subclasses. 
Considering this fact, the model did surprisingly veil. All class 9000 
(me tail structures) vere correctly classified, a 100JS performance. 

6.4 CONCLUSIONS 

Since tlie results presented here are preliminary , no optimization or 
detail analysis of the model or variants has been accomplished . Hov— 

ever , it .1 ; believed that the evidence presented in this paper indicates 
that map to map scintillations for the type SAR studied do contain potential- 
ly valuabli object and feature classifier information . 
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